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Abstract 

We consider the open spin-s XXZ quantum spin chain with nondiagonal boundary 
terms. By exploiting certain functional relations at roots of unity, we propose the 
Bethe ansatz solution for the transfer matrix eigenvalues for cases where atmost two 
of the boundary parameters are set to be arbitrary and the bulk anisotropy parameter 

has values r\ = ^ , ^ , We present numerical evidence to demonstrate completeness 

of the Bethe ansatz solutions derived for s = 1/2 and s = 1. 
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1 Introduction 



There have been significant focus of effort in solving integrable quantum spin chains for 
many years. In particular, integrable quantum spin chains with boundaries (integrable open 
quantum spin chains) have attracted much interest over the years. As a result, models such 
as the open XXX and XXZ quantum spin chains have been subjected to intensive studies 
due to their growing applications in various fields of physics, e.g. statistical mechanics, 
string theory and condensed matter physics. Despite numerous success in the past p]-[6] 
(also refer to [T]-[ll] and references therein, for other related work on the subject.), there 
still remain unsolved problems in this area. Bethe ansatz (in its conventional form) for the 
most general case of the open XXZ quantum spin chain (even for the spin-1/2 case) with 
arbitrary nondiagonal boundary terms and generic bulk anisotropy parameter is yet to be 
found. In [14], Galleas found an interesting solution analogous to Bethe ansatz equations 
for the spin-1/2 case. This solution, written in terms of certain functional relations are 
expressed in terms of roots of the transfer matrix. Much progress have been made on the 
topic up to this point. In a series of publication, Bethe ansatz solutions have been derived 
for open spin-1/2 XXZ quantum spin chain where the boundary parameters obey certain 
constraint. Readers are refered to [15] -[19] for related work on the subject. Apart from this 
constraint, two sets of Bethe ansatz equations are needed there to obtain all 2 N eigenvalues, 
where N is the number of sites. A special case of the above solution was generalized to 
open XXZ quantum spin chain with alternating spins by Doikou [20] using the functional 
relation approach, proposed by Nepomechie in [16] to solve the spin-1/2 case (which indeed 
the method used in this paper). In [21], related work was carried out using the method in 
|15j . Recently in [22], Frappat et al. further generalized the spin-1/2 XXZ Bethe ansatz 
solution (for boundary parameters obeying the constraint) to the spin-s case by utilizing 
an approach based on Q-operator and T-Q equation, which was developed earlier for the 
spin-1/2 XXZ chain in [TH] and subsequently applied to the spin-1/2 XYZ chain in [23] . As 
in the spin-1/2 case, two sets of Bethe ansatz equations are also needed there to produce all 
(2s + 1)^ eigenvalues, where again N represents the number of sites. 

In this paper, we present Bethe ansatz solutions for open spin-s XXZ quantum spin chain 
without such a constraint among the boundary parameters. We follow similar approach as 
given in [161 El 121] that was used to solve the s = 1/2 case. It is based on fusion [U [25], 26], 
the truncation of the fusion hierarchy at roots of unity [27] and the Bazhanov-Reshetikin 
[28] solution of the RSOS models. As in [21], there are atmost two arbitrary boundary 
parameters. The rest of the parameters are fixed to some values. The approach we use, 
which is based on functional relations obeyed by transfer matrix at roots of unity [16] yields 
Bethe ansatz solution which gives completely all the (2s + 1)^ eigenvalues. One limitation 
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of the solution is that it is valid only at roots of unity, namely when the bulk anisotropy 
parameter has values 77 = r^j- In this paper, we consider only even values of p. Lack of single 
set of Bethe ansatz equations that yield complete eigenvalues for the model considered here, 
namely where the boundary parameters are arbitrary (even at most two) has motivated us to 
study this problem. Moreover, we note that the relation of s — 1 case to the supersymmetric 
sine-Gordon (SSG) model [29J (here the boundary version [SU] EI]), has also been part of 
our motivation for considering the problem. 

The outline of the paper is as follows: In Sec. 2, we review the construction of the 
so-called fused R [25], [32] [33] El] and jH [26] matrices from the corresponding spin- 1/2 
matrices. For some original work on spin- 1/2 K T matrices, refer to [35] [36] • Construction 
of commuting transfer matrices from these fused matrices (using Sklyanin's work [3J, which 
in turn relies on Cherednik's previous results [M]), together with some of their properties 
are reviewed. Fusion hierachy and functional relations obeyed by transfer matrices are also 
reviewed. In Sec. 3, we present the Bethe ansatz solutions for cases with atmost two 
arbitrary boundary parameters at roots of unity, e.g. 77 = y , y , . . ., by exploiting the 
reviewed functional relations obeyed by the transfer matrices. Further, we present numerical 
results in Sec. 4 to illustrate the completeness of our solution, using s = 1/2 and s = 1 
as examples, where the Bethe roots and energy eigenvalues derived from the Bethe ansatz 
equations (for some values of p and N) are given. We remark that these energy eigenvalues 
coincide with the ones obtained from direct diagonalization of the Hamiltonians. Finally, we 
conclude the paper with discussion of the results and potential future works in Sec. 5. 



2 Transfer matrices, fusion hierachy and functional re- 
lations at roots of unity 

In this section, in order to make the paper relatively self-contained, we review some crucial 
concepts on the construction of commuting transfer matrices for iV-site open spin-s XXZ 
quantum spin chain. Materials reviewed here on fused R, K T and higher spin transfer matri- 
ces are borrowed from [22], as presented there. As constructed in [3J, the commuting transfer 
matrix for s = 1/2, which we denote (following notations adopted in |22j ) by t^'^(w), whose 
auxiliary space as well as each of its N quantum spaces are two-dimensional, one can similarly 
construct a transfer matrix t^' s \u) whose auxiliary space is spin-j ((2j + l)-dimensional) and 
each of its iV quantum spaces are spin-s ((2s + l)-dimensional), for any j, s G {|, 1, |, . . .} 
using the so-called fused R [25] E21 [33] El] and K T jH|26] matrices. As for the spin-1/2 case, 
these R and K T matrices serve as building blocks in the construction of the commuting 
transfer matrices for higher spins. We list them below along with some of their properties. 
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The fused- R matrices can be constructed as given below, 



ftw = p w p w n n (« + (* + 1 - 3 - s - !)»?) p w p < 



{6} 



(2.1) 



fe=i i=i 



where {a} = {ai, . . . , a2j}, W = {&i> • • • , &2s}> and P^, is the symmetric projector given by 



1 23 ( k 

T.> = OT II(E* 

V fe=i \z=i 



(2.2) 



"P is the permutation operator, with V ak ,a k = 1; Similar definition also holds for P ( 
(u) is given by 



( sh(u + 77) \ 

shw sh rj 

shrj shu 

\ sh{u + r]) ) 



(2.3) 



where 77 is the bulk anisotropy parameter. Note that the fundamental R matrix satisfies the 
following unitarity relation 



R^^\u)R { ^\-u) 



-£(«)! , = sh(w + rf) sh(u — rf) . 



(2.4) 



The R matrices in the product (12.11) are ordered in the order of increasing k and /. The fused 
R matrices satisfy the Yang-Baxter equations 



(2.5) 



- «) = 4*4}^) ftw *&V u - v ^ ■ 



Having defined fused-R matrices, one can analogously construct fused K matrices [U [26] 

rfc-i 



>} 



2j 



k=l 



l[RiUr(2u+(k + l-2j-l)r 1 ) 

1=1 



x ^ ) (n + (A;-j-i)r / )}p { + a} , 

where K~^\u) is the 2x2 matrix whose components are given by 

K\x ( u ) = 2 (sh a_ ch (3— ch u 4- ch ct_ sh /3_ sh w 
K^iu) = 2 (sh a_ ch /?_ ch u — ch «_ sh /3_ sh w 
lff 2 (u) = e e - sh 2u , K^(u) = e~ 9 - sh 2u , 



(2.6) 



(2.7) 
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where a_ , /3_ ,0_ are the boundary parameters. The products of braces {. . .} in (12.61) are 
ordered in the order of increasing k. The fused K~ matrices satisfy the boundary Yang- 
Baxter equations [38] 



{a} V / {a}{b}V 



'{6} V"; "{«}{&} 

The fused matrices are given by 



1 



where the normalization factor is. 



^a} } (-«-^) 



2J-1 I 

f u \u) = n n[-^ 2u + v + fc + 1 - 2 ^»] 

Z=l fe=l 

Using the above results, one can construct the transfer matrix t^' s \u), 

*<">(«) = tr {a} tf^>(«) f { ^(«) , 

where the monodromy matrices are given by products of iV fused R matrices, 



{21 



(2.9) 



(2.10) 



(2.11) 



R 



{a},{bim } 
l {a},{6W}^ 



V)... i2 



0» 

{a},{6W} 



:«) 



l {a},{bl N l} [ 



(2.12) 



These transfer matrices commute for different values of spectral parameter for any j,j' e 
{|, 1, |, . . .} and any s G {|, 1, §, . . .}, 







(2.13) 



Furthermore, they also obey the fusion hierarchy [U [26], 22]^. 



t (i-i^)( M _ J7] ) t ^ s \u) = t^'\u - (j - ^)v) + S {s) {u - V ) t^'\u - (j + -)r}) , (2.14) 



j = 1, |, . . ., where t^°' s ^ = 1, and 5^(u) is given by 



(»)/ 



(2.15) 



where 



'2s- 1 



fc=0 



2JV 



sh(2w) sh(2« + Arj) 



sh(2w + 77) sh(2w + 3t?) 

8[ (u) = 2 4 sh(u + a_ + 77) sh(w — a_ + 77) ch(u + /?_ + 77) ch(tt — /3_ + 77) 

x sh(-u + a + + 77) sh(w — a + + 77) ch(w + (3 + + 77) ch(w — (3 + + rf) . (2.16) 



^ee the appendix in |22j for more details on the fusion hierachy. 



Note that the 5^ (u) in [22] differs to the one given here merely by a shift in 77. 

Next, we list few important properties of the rescaled "fundamental" transfer matrix 
t^' s \u) (defined below), which are useful in determining its eigenvalues. Following the 
definition of t^ ,s \u) as in [22], we have 



t^ s \u)= ^ — tG*\u), (2.17) 

where 



2s-l 

g (h»)(u) = Jl sh(u + (s-k + ~)rj) (2.18) 

k=l 

This transfer matrix has following useful properties: 

i { ^ s \u + m) = t^' s) (u) (in - periodicity) (2.19) 

t { ^ s \-u - 77) = t { ^' s) (u) (crossing) (2.20) 
f(i s )(0) = -2 3 sh 2iV ((s + ^)r/)chr/sha_ch/5_sha+ch/5+I (initial condition) (2.21) 



= 2 3 sh 2N u 

r;=0 



— sh a_ ch P_ sh a + ch (3 + ch 2 u 
+ ch a_ sh P_ ch a + sh (3 + sh 2 u 

- ch(#_ -# + )sh 2 tich 2 w I (semi-classical limit) (2.22) 



f*>'\u) -l-e^+^+^+^ch^-^)! for u — > +00 

(asymptotic behavior) (2.23) 

where I is the identity matrix. 

Due to the commutativity property (12.131) . the corresponding simultaneous eigenvectors 
are independent of the spectral parameter. Hence, (12.191) - (12.231) hold for the corresponding 
eigenvalues as well. In addition to the above mentioned properties, for bulk anisotropy values 
77 = ^e_, with p = 1,2,..., the "fundamental" transfer matrix, t^' s \u) (and hence each of 
the corresponding eigenvalues, A^' s \u)) obeys functional relations of order p + 1 

t {1 i' s \u)t^> s \u + 77) . . . t&'\u + pq) 
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- 5^ s \u - r])t^' s \u + r])t^' s \u + 2rj)... t^' s \u + (p - 1)77) 

- 5 {s) {u)t^' s \u + 2 V )t^' s \u + 377) . . . t^' s \u + prj) 

- 5 {s \u + ri)t { ^ s \u)t^> s) (u + 3r])t^' s) (u + A-q) . . . t& a) (u + prj) 

- 5 {s \u + 2r l )t^ s \u)t^ s \u + rj)t^' s \u + 477) . . . S' s \u + prj) - . . . 

- <5 w (w + (p - l)?7)^' s) (u)^' s) (u + 77) . . . ^' s) (u + (p - 2)77) 

+ ... = /(«). (2.24) 
For example, for p = 2 and p = 4, the functional relations are 

*(3> s )(«)*(2' s >(« + r])t ( h s \u + 27/) - 5 (s) (w - rj)t& a) (u + 77) - (jW(u)^3> 8 )(u + 277) 
-5 (s \u + 77)^^) = f(u) . (2.25) 

and 

t(i s )( M )t(i s )( M + r])S' s \u + 2ry)7i ( ^ s) (n + 3r])t^' s) (u + 4t/) 



-<5<« 



} (w + r])5 (s \u - 2rj)t^ s \u) + 5 (s) (w)5 (s) (w + 277)^2.-) ( u + 477) 
\u + r])5 (s) {u - rj)t^' s \u + 3r]) - 5 {s) (u + rj)t^ s) {u)t^' s) {u + 3r])t^' s \u + Arj) 
} (m)<5 (s) (m - 2r])t^' s) (u + 277) - 5 {s \u)t^' s) (u + 2r])t^ s \u + 3r])S' s \u + 477) 
} (« - rj)5 {s \u + 277)t ( 5- s )( M + 77) - 5 (s) (w + 27])t^' s) (u)t^' s) (u + 7])t^' s) (u + 4 V ) 
\u - 2r])t ( ^ s \u)t^' s \u + r])t^' s \u + 277) 

\u - r])t^' s) (u + rj)t^' s \u + 2rj)t^' s) (u + 377) = f(u) . (2.26) 



respectively. The scalar function f(u) (which can be expressed as f(u) = f (u)fi(u)) is given 
in terms of the boundary parameters a T , f3 T , 6 T (for even p) by 



fo(u) 



-l) N+1 2- 4s P N sh 4sN ((p + l)u) , 
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;-i) 7V+1 2- 4 ^ 7V ch 4sAr ((p+ i)u) , 

S = l,2,3,... 



(2.27) 



and 



fi(u) = (-lf +1 2 3 ^( 

sh ((p + l)a_) ch ((p + sh ((p + l)a+) ch ((p + 1)/?+) ch 2 ((p + l)u) 

- ch ((p + l)a_) sh ((p + ch ((p + l)a+) sh ((p + sh 2 ((p + l)u) 

- (-l) iV ch((p+l)(^-^ + ))sh 2 ((p+l) M )ch 2 ((p+l) M )). (2.28) 



sh ((p + l)a_) ch ((p + sh ((p + l)a+) ch ((p + 1)/?+) ch 2 ((p + 

- ch ((p + l)a_) sh ((p + ch ((p + l)a+) sh ((p + 1)0+) sh 2 ((p + 1)«) 

- ch((p+l)(^-^ + ))sh 2 ((p+l) M )ch 2 ((p+l) M )). ( 



(2.29) 



for s = 1 , 2 , 3 



. . ., Note that f(u) satisfies 



f(u + r,) = f(u), 



f(-u) = f(u). 



(2.30) 



and 



p 



(2.31) 



i=o 



where 6q (u) is given by (12.161) . 

3 Bethe ansatz 

In this section, we give main results of this paper. We derive Bethe ansatz equations for 
various cases where atmost two of the boundary parameters {a_, a+, /?+} are arbitrary 
by adopting the steps given in [23]. By considering atmost two boundary parameters, we find 
certain factors in the calculation become perfect squares. This facilitate the computations 
that follow. More on this is explained below. 

3.1 a + , a- arbitrary 

Here, we take both a_ and a + to be arbitrary while setting (3± = rj, = 9 + = 6, where 9 
is arbitrary. In order to obtain Bethe ansatz equations for the transfer matrix eigenvalues 
Ail') (u), we shall recast the functional relations (12. 24ft as the condition that the determinant 
of a certain matrix vanishes (following [28J). We find that the functional relations (I2.24p for 
the transfer matrix eigenvalues can be written as 



det M = , 



(3.1) 
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where M. is given by the (p + 1) x (p + 1 ) matrix 



M 



( A^' s \u) -h(u) 

—h{—u) A^' s \u + prj) —h{u+prj) ... 








-h(—u + pr]) ^ 




i —h{u + p 2 rj) 











—h(—u—p(p—l)r)) A.^' s \u + p 2 T]) j 



(3.2) 



(whose successive rows are obtained by simultaneously shifting u t— > u + pr] and cyclically 
permuting the columns to the right) provided that there exists a function h(u) with the 
following properties 



h(u + 2m) = h(u + 2(p + l)rj) = h(u) , 

h(u + (p + 2)77) h(-u -{p + 2)77) = (u) 
p p 

Y[h(u + 2 JV )+l[h(-u-2j V ) = f(u). 

3=0 j=o 



(3.3) 
(3.4) 

(3.5) 



From (I3.3l) -( l3~5l) . we see that the problem of finding h(u) then reduces to solving the following 
quadratic equation in z(u), 



z\u) — z\u 



)/(«) + fp W (u + (2j - 1)^) = 

j=0 



(3.6) 



where 



p 

z{u) = Y[ h{u + 2j V ) . 

j=0 



(3.7) 



For the cases considered here and in subsequent sections, the discriminants of the corre- 
sponding quadratic equations are perfect squares, and the factorizations such as (13.71) can 
be readily carried out. However, when all boundary parameters are arbitrary, the discrimi- 
nant is no longer a perfect square; and factoring the result becomes a formidable challenge. 
Solving the quadratic equation (13.61) for z(u), making use of the explicit expressions (12.161) 
and (I2.27p - (l2.29p for 5^ s \u) and /(«), respectively, we obtain the following for h(u), 



with 



ho(u) = (-iy sN A 



h{u) = h (u)hi(u) , 



2s-l 

i J sh(u + (s - k + -)r]) 

.k=0 



(3.8) 



2iV 



sh(2w + 277) 
sh(2w + 77) 



(3.9) 
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and 



hi(u) = < 



ch 2 (u - n) sh(u - a.) sh(u + ( _ 1)j v a+) ^(j(^-H»)^(^( M+ (-i)^ a++ ,)) 
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2 ' 2 ' 2 ' - - - 



ch 2 (u — rj) sh(u — a_) sh(u + a + 



ch(§(-u-a_-r?)) ch(±(«+(-l) Jv a + -» ? )) ' 

ch(i(u+a_+7y)) ch(i(«-a++r?)) 
ch(|(u— a_— rj)) ch(i(u+a:+— »?)) ' 



(3.10) 



S = l,2,3,... 



Further, the structure of the matrix M. (13.21) suggests that its null eigenvector has the form 
(Q(u) ,Q(u + prj) , . . . ,Q(u + p 2 r])), where Q(u) has the periodicity property 



Q(u + 2m) = Q(u) . 
It suggests that the transfer matrix eigenvalues are given by 

A ( 2' s) (u) = h{u) r + h(-u+pq)- 



Q(u) 



Q{u) 



(3.11) 



(3.12) 



which is of the Baxter's TQ relation form. Noting that the functions h(u) and h{—u+prf) (see 
fl3jP-( [3Jl)D ) have the factor g^ s \u) 2N (see (EE])) in common (since g { ^ s \v) = g&'\-u + 
prj)), we can rewrite f 1 3 . X 2 1) in terms of the eigenvalues of t^ ,s \u) (see (I2.17P ) as 



A<**>(«) = + H-u+pr,)^-^ 



Q(u) 



Q(u) 



where 



(3.13) 



with 



and 



h{u) = ho(u)hi(u) 



h Q (u) = (-l) 2sN A sh 2N (u + (s + ^77) 



1. ,sh(2u + 2r/) 



2' " sh(2u + r]) 



Q(u) =J[sh Q(«-t*)) shQ( 



■u + Uj — prj) 



with the periodicity (13. lip as well as the crossing property 

Q(-u + prj) = Q(u) . 

where 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



M = 2sN + 2p + 1 



(3.18) 
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which is confirmed numerically for small values of N and p. We stress here that the h(u) 
given above is not the only solution. It is obtained largely by trial and error, verifying 
numerically for small values of TV that the eigenvalues can indeed be expressed as (13.131) 
with Q(u) y s of the form given by (13.161) . We also remark that (13 . 18[) is consistent with the 
asymptotic behavior (I2.23p . Making use of the analyticity of A(2' s \u), we have the following 
for the Bethe ansatz equations, 

h{uj) Q{uj-pr]) 



h(-Uj+pr]) Q{uj+pri) 



M. 



(3.19) 



3.2 (3 + , f3- arbitrary 

In the following, we set /?_ and f3 + arbitrary while setting a± = rj, 6L = 9 + = 9. As before, 
we write the functional relations ( 12.241) for the transfer matrix eigenvalues in the form of 
( 13. ip . where for this case, the matrix M. is given by 

M = 

( A^' s \u) -h{u) ... 

— h(— u — (p + 1)77) A^' s \u + prj) —h(u+prj) ... 



—h(—u — 77) ^ 




y —h{u + p 2 rf) 











-h(-u - (p 2 + 1)77) A ( i s) (w + p 2 i]) J 



(3.20) 



if h(u) satisfies 



h(u + 2m) = h (u + 2(p + 1)77) = h(u) , 

h(u + {p + 2)77) h(-u - 77) = 5 {s) (u) 
p p 

Y[h(u + 2jr ] ) + l[h(-u-(2j + l)r ] ) = f(u). 

j=Q j=0 



(3.21) 
(3.22) 



(3.23) 



Proceeding in a similar way to the previous case and setting h(u) = ho(u)hi(u) we find 



ho(u) = (-l) 2sw 4 



2s-l . 

il sh(u + (s - k + -)v) 



k=0 



2N 



sh(2w + 2r]) 
sh(2w + 77) 



(3.24) 



and 



hAu) 



sh(u — rj) sh(w + rj){ch.u — % sh /3_)(chw + (— l) N i sh (3+) , 

e = I 3 5 
^ 2 , 2 , 2 , • • • 

sh(w — 77) sh(w + ?7)(chM + i sh /?_)(chw — i sh j3 + )} , 



(3.25) 



S = l,2,3,... 
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The transfer matrix eigenvalues are now given by 

A<M(„) = Ku)9t±^ + - , (3.26) 

Q{u) Q{u) 

As before, due to the common factor g(h s )( u ) 2N (see 12,181) , and using the crossing symmetry 
gr(|' s )(u) = ±g(i' s )(— u — 77), we conclude that the eigenvalues of t^ ,s \u) are given by 

A<i.->(„) = M«)^±M + ft(-u - , (3.27) 

Q(u) 



where 



and 



h(u) = /ioH^i(m) (3.28) 



»„(„) = (_l)-4 dl »(« + (. + I),)^±M ,3.29) 



The ansatz for Q(w) is given by 

Q(u)=J[sh(^( U -u j ^sh(^(u + u,+ V )^ , (3.30) 

which satisfies Q{u + 2«7r) = Q(u) and Q(—u — rj) — Q{u); and 

M = 2siV + p. (3.31) 
Moreover, the Bethe ansatz equations for the zeros Uj take the form 

h(uj) _ Q(uj-pq) 



h(-Uj-rj) Qiuj+pv)' 



j = l,...,M. (3.32) 



where we find the number of Bethe roots (13.311) is consistent with the asymptotic behaviour 
3.3 One arbitrary (3 and one arbitrary a 

Finally, we consider combinations where the arbitrary parameters consist of one of the /3's 
and any one of the a's. To keep the expressions general, we drop the subscripts ± from the 
boundary parameters, a±, (3±. The remaining boundary parameters are fixed, e.g., 
arbitrary, f}_ = 77, a + = y or other similar combinations. Also, as in previous cases, we let 
6L = 9 + = 9. The matrix M. is identical in form as in (13. 2D . 
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We once again find h(u) = h (u)hi(u), with the same Hq{u) as for the earlier cases. For 
hi(u), we take the following, 



hAu) 



ch u ch(u — 77) (sh u + (—l) N i ch 8) sh(u — a) ch {^ u+a+v ^) 

' ch( ii(u— a— rf)j 



o _ I 3 5 
15 2 ' 2 ' 2 ' 



ch u ch(w — 7]) (sh u + idid) sh(w — a) -44 



ch(i(u+a+?7)) 



(3.33) 



di(|(u-a-7?)) 



S = l,2,3, 



The above /i(u) satisfies (I3.3l) - (l3.5p . The eigenvalues of the transfer matrix and Bethe ansatz 
equations are given by A3. 12j) . (13.131) . (13. 16|) and (13.191) . with 

M = 2Ns+p (3.34) 

which again is consistent with (12.231) . We note that for s = 1/2, our solutions for all the 
above cases coincide with the corresponding solutions found in 



4 Energy eigenvalues and Bethe roots 

In this section, we illustrate the completeness of the Bethe ansatz solutions derived in Sec. 
3. We provide numerical evidence for cases s — 1/2 and s = 1, namely the complete energy 
levels together with the Bethe roots used in the computation (see Tables 1 and 2). 

4.1 s = 1/2 case 

The Hamiltonian for the open spin-1/2 XXZ quantum spin chain is given by [351 EH] 



n = 2Ete + *»+i +ch? ' (7 «<+i) 

n=l 

+ — sli ?7 T coth a_ tanh f3-o{ + csch a_ sech 8-{ ch Q^a\ + % sh 9-crf) 

— coth a + tanh [3 + a z N + csch a + sech 8 + ( ch^+a^ + i sh 6 + a v N ) , (4.1) 

where a x ,a y ,a z are the standard Pauli matrices, rj is the bulk anisotropy parameter, 
a ± 5 P± 5 Q± are arbitrary boundary parameters, and N is the number of spins. 

We compute the energy eigenvalues of (14.11) (from Bethe ansatz) for a particular case 
derived in Sec. 3. For the purpose of illustration, it is sufficient to consider the case where 
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the two arbitrary boundary parameters are a_ and /3_. The steps here can be repeated for 
any other desired combinations of boundary parameters^ The Hamiltonian ( 14 .11) is related 
to the first derivative of the transfer matrix, t^>|)(^)[3] 



du 



u=0 



4^1, (4.2) 



where 



16 sh a_ ch (3- sh a + ch f3 + sh 27V 1 rj ch 77 



(i) sh 2 77 + iV ch 2 n 

2 ch?7 

and I is the identity matrix. Moreover, f!4.2[) implies that the energy eigenvalues are given 
by 

E = + c { }\ (4.4) 

du u=o 

Hence, using the results (13 . 13[) - (13 . 16[) and (l3.33ll R one arrives at the following result for the 
energy eigenvalues in terms of Bethe roots {uj}, 

rn 1 1 1 V 1 1 at 1 1 ch2?7 

hi = -sn?7cn— > — —, — -r- — H — J\ can 

2 ' 2^sh(f« i )ch(f(w,- + 77)) 2 ' 2 ch V 

1 01— + 11 

— — sh^(cotha;_ + % sech /5_ — tanh( — - — )) . (4.5) 

where M = N + p (see ( l3~34j) l 

In Table 1, we tabulate the energy eigenvalues computed using (14.51) for = 4 together 
with the Bethe roots (These roots are obtained using a method developed by McCoy and 
his collaborators [3H] which is also explained in [IT].). This numerical result illustrates the 
completeness of Bethe ansatz equations derived in Sec. 3. We have verified that the energies 
given in Table 1 coincide with those obtained from direct diagonalization of (14.11) . 



4.2 s = 1 case 



In this section, we repeat the analysis for s = 1. We shall consider the case investigated in 
Sec. 3.2, namely the case with arbitrary (3 + . The integrable Hamiltonian for the open 
spin-1 XXZ quantum spin chain is given by (adopting notations used in [22J) 

N-l 

H=J2 H n,n+l + H b . (4.6) 
n=l 

2 Note that for s = 1/2, t^'i\u) = t^^\u) 

3 The function h(u) used here coincides with the one found in [23] , 
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H n;H+ i represents the bulk terms. Explicitly, these terms are given by [3D] , 
H n , n+1 = a n - (a n f + 2 sh 2 77 [< + {S z n f + (S z n+1 ) 2 - «) 2 ] 



where 



4sh 2 (|) + 



&n — S n ■ S n+ \ , a n — S n S n+1 + S^S^ +t , a n — S n S n+l 



(4.7) 



(4i 



and S" are the sw(2) spin-1 generators, if;, represents the boundary terms which have the 
following form (see e.g., 



(4.9) 



H b = cuiSl) 2 + a 2 S{ + a 3 {S+) 2 + a 4 {S^) 2 + a 5 S+ S z t + a 6 S* Sf 
+ a 7 Sl Sf + a 8 Si Si + (aj <h> bj and 1 <-> N) , 

where = S x ±iS y . The coefficients {a{\ of the boundary terms at site 1 are functions of 
the boundary parameters (pt-,/3-, 6L) and the bulk anisotropy parameter 77. They are given 
by, 

1 

T 
1 

I' 

1 



a i 



(i-2 



-Oq (ch 2ct_ — ch 2/3_ + ch 77) sh 2rj sh 77 . 
-cio sh 2a_ sh 2(3- sh 277 , 



«3 

f/7 
«8 



--aoe 2e_ sh277sh77 , 
8 

— a e~ 26 ~ sh 2?7 sh 77 , 
8 

aoe 61- ( ch /3_ sh a_ ch — + ch a_ sh (3_ sh — J sh 77 ch 5 77 , 
aoe -6- ^ch /3_ sh a_ ch — + ch a_ sh /3_ sh — J sh 77 ch 5 77 

/ 7^ 77 \ 

— a^e 6 " ych. /3_ sh a_ ch — — ch a_ sh /3_ sh — J sh 77 ch 
— aoe - ' 



77 



where 



a 



77 77 s 

ch /3_ sh a_ ch ch a_ sh /3_ sh - ) sh 77 chs 77 . 

2 2 / 



sh(a_ - |) sh(a_ + |) ch(/3_ - |) ch(/L + | 



(4.10) 



(4.11) 



Similarly, the coefficients of the boundary terms at site iV which are functions of the 
boundary parameters (a+, /3+, 9 + ) and 77, are given by the following correspondence, 



bj = a, 
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(4.12) 



To derive the energy formula similar to (14.51) for s = 1 case, we once again begin by 
expressing the spin-1 Hamiltonian in terms of the first derivative of spin -1 transfer matrix, 
namely t^'^(u). One can construct t^(u) from t^' 1] (u) by using the fusion hierarchy 
formula (jHHD , 



t^\u)=t^\u-l)t^\u + l)-^\u~l) 



(4.13) 



where 8^ l '{u) is given by (I2.15l) - (l2.16p with s = 1. Following |22j . we work with rescaled 
transfer matrix given by 



t 



(i,i) gt. 



sh(2u) sh(2u + 2 V ) f(1>1) 



(4.14) 



[shw sh(w + ?7)] 27V 

where t*- 1 ' 1 ) 9t (u) is the transfer matrix constructed from "gauge" -transformed R^^u) and 
K T ^(u) matrices 0. Note: The rescaled transfer matrix does not vanish at u = 0. The 
Hamiltonian H (14. 6ft . according to [3], is related to the first derivative of ft 1 ' 1 ' 9t (u), 



du 



u=0 



(4.15) 



which in turn implies that the energy eigenvalues in terms of transfer matrix eigenvalues 
A.(i>i) 9 *(w), are given by 



du 



u=0 



(4.16) 



where 



and 



c (1) 



ch^|l6[sh2^sh?7] 27V sh3r/sh(a_ — — ) sh(a_ + — ) ch(/3_ 



x sh(a+ _ |) sh(a+ + 1) ch(/3+ _ |) ch (/3 + + 1; 



-^6 ch - (N - 1) (4 + ch 2rj) + 2iV ch 2 rj 



(4.17) 



shr? 
~2d 



- 2ch2a + (chr?(3 + 7ch2?7 + ch4r/) + ch2/3+(4 + 5ch2r/ + 2ch4?7)) 

+ 2chr/^ch2/3+(3 + 7ch2r/ + ch4r/) + chr/(5 + 3 ch2r? + 3 ch4r?) J } 

sli 2?7 r / 

f 1 ch2/?+(2 + 4ch?7ch377) + chr](5ch2r] + ch47?) - 2 ch 2a+ ( 1 + ch2?7 



2d 



+ 



ch 2/5+ (ch 77 + 2 ch 377) + ch 4r]) j 



(4.18) 



One reason for such a transformation is to bring these matrices to a more symmetric form. For a detailed 
discussion on this, refer to Sec. 4 of 22J 
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where 



and 



b = 2( - ch2p_ - ch 3 77 + ch2a_(l + ch2/?_ chr])) (4.19) 



d = -4 sh 377 sh(a + + |) sh(a + - |) ch(/3 + + |) ch(/3 + - |) (4.20) 

Furthermore, using the fact that A (1>1) 9t {u) = A^'^ii) and (pT27jl . (Q3j) . (14141 . we obtain 
the energy in terms of Bethe roots {uj}, 

£ = --shr ? sh2r ? ^ — + - _ , ^ (A'(0) + i?'(0)) 

2 ^ch(i(^ + ^))ch(i( Mj -2)) 2^2)/ l (_2) 

+ c^C'^+cW. (4.21) 

where 

A( W ) = fc( u +|)X(ti-|) 

B{u) = -h{u+\)k-u-\) 

ch(w + 2) -ish/3_^ ^ch(w+ 2) +ish/? 4 



(( '° " \ch( M +2)+zsh/3_Kch( M + f)-^ s h/?J /,>! ''" 



= 4 sh 2iV (w + y ) sh(2w + 2rj) sh(w + 77) sh(w - 77) 
x (chw + zsh/3_)(chu-zsh/3+) (4.22) 

Also, M = 2N + p (see (13T3T]) ). 

We tabulate the energies computed using (14.211) for AT = 3 with the Bethe roots (which 
are obtained using similar method as for the s = 1/2 case above) in Table 2. These numerical 
results once again illustrate the completeness of Bethe ansatz equations derived in Sec. 3. 
We have verified that the energies given in Table 2 coincide with those obtained from direct 
diagonalization of (14.61) . One can proceed to repeat the analysis for higher spin values, 
namely s > 1. However, due to tedious computations, we avoid from pursuing it here. 



5 Discussion 

We have determined Bethe ansatz solutions of the open spin-s XXZ quantum spin chain 
for cases with nondiagonal boundary terms (I3.13p - (l3.19p and (I3.27I) - (I3.32I) . by following the 
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method used earlier in [16, 24J to solve the spin-1/2 case. This method relies on functional 
relations (12.24ft that the "fundamental" transfer matrices, t^ ,s \u) obey at roots of unity. 

However, these solutions hold only for 77 = y , y , Unlike Bethe ansatz solutions found in 

earlier works on the open spin-s XXZ chain with nondiagonal boundary terms, we emphasize 
that Bethe ansatz solutions found here hold for arbitrary values of boundary parameters 
(atmost two). We have checked these solutions for chains of length up to iV = 4, and have 
verified that indeed they give the complete set of (2s + 1)^ eigenvalues. Moreover, we also 
presented numerical evidence for the completeness of the Bethe ansatz solutions found (using 
s = 1/2 and s = 1 as examples) in Tables 1 and 2. Perhaps the completeness of the Bethe 
ansatz equations for spin-s can readily be concluded from completeness of the corresponding 
Bethe ansatz equations for spin-1/2 case and the fusion hierarchy (12. 14ft which is used in the 
construction of higher spin-s transfer matrices. 

There remain many problems worth investigating. As mentioned in the Introduction, 
due to the relation of s = 1 case to supersymmetric sine-Gordon (SSG) model, one can carry 
out similar analysis as in [31], but now for spin-1 chain with nondiagonal boundary terms. 
One could also try to extend the solutions presented here to cases with multiple Q(it)'s as in 
[42] . Also, to our knowledge, conventional form of Bethe ansatz solution for the open XXZ 
quantum spin chain, where all six boundary parameters are arbitrary with generic values 
of bulk anisotropy parameter 7], has not been found. We remark that through a series of 
important work on the spectrum of XXZ spin chain based on representation theory of the 
q-Onsager algebra [TT], Baseilhac and Koizumi argue that obtaining such a conventional 
Bethe ansatz solution for the most general case is unlikely. It would be interesting to explore 
their results further and compare their approach with the Bethe ansatz approach. 
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23 



E 



Bethe roots «, 



-3.19769 



-2.42188 



-1.87006 



-1.30053 



-0.874711 



-0.674656 



-0.203476 



0.222018 + 2.91719 i, 0.0900395 + 2.91719 i, - 2.6018 i, 
1.01834 - 1.7952 i, 2.15279 i, 0.267003 - 1.7952 i, 
1.00769 i, 1.0165 + 1.3464 i, 0.0900395 - 0.224399 i, 
0.222018- 0.224394 i 
0.324807 - 3.13487 i, 0.319576 + 2.68662 i, 0.0958764 + 2.91717 i, 
-2.60637 , 2.15279 i, 0.356519 - 1.7952 i, 
1.09238 i, 0.0958764 - 0.224378 i, 0.324807-0.455523 i, 
0.319576 + 0.00617674 i 
0.530712 + 2.91722 i, 0.0853747 + 2.91719 i, -2.60166 i, 

0.946517 - 1.7952 i, 2.15279 i, 0.25634 - 1.7952 i, 
1.00255 i, 0.943646 + 1.3464 i, 0.0853747 - 0.224399 i, 
0.530712 - 0.224428 i 
0.0805934 + 2.91719 i, 1.4454 - 1.7952 i, -2.6016 i, 
0.607877 - 1.7952 i, 0.212672 - 1.7952 i, 0.592602 + 1.3464 i, 
0.992082 i, 0.54 i, 1.44531 + 1.3464 i, 
0.0805934 - 0.224399 i 
0.284541 + 2.91781 i, 0.251638 - 3.13998 i, 0.245572 + 2.69145 i, 
2.15279 i, 0.359151 - 1.7952 i, 1.59561 i, 
-0.983447 i, 0.251638 - 0.450414 i, 0.245572 + 0.00134169 i, 
0.284541 - 0.225016 i 
0.518223 + 2.91722 i, 0.199566 + 2.91719 i, 0.940199 - 1.7952 i, 
2.15279 i, 0.255178 - 1.7952 i, 1.00201 i, 
-0.988749 i, 0.937209 + 1.3464 i, 0.199566 - 0.2244 i, 
0.518223 - 0.224431 i 
0.182373 + 2.91719 i, 1.43967 - 1.7952 i, 0.604391 - 1.7952 i, 
0.211794 - 1.7952 i, 0.588802 + 1.3464 i, 0.991943 i, -0.988797 i, 
0.54 i, 1.43958 + 1.3464 i, 0.182373 - 0.224399 i 



Table 1: The 16 energies and corresponding Bethe roots given by A^'^'(n 



N = 4,s 
ztt/2 , /3+ = 



= 1/2, p = Q,v 
V,9 + = 



in/7, a_ = 0.54i = 0.2, 9. 



\u) for 
, a + = 
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E (continued) 



Bethe roots Uj (continued) 



0.343441 



0.761262 



0.846541 



0.883622 



1.00689 



1.20648 



1.50502 



1.82374 



2.16601 



0.149446 - 3.14159 i, 0.149313 + 2.69279 i, 1.05439 - 1.7952 i, 
2.15279 i, 0.273366 - 1.7952 i, 1.01147 i, 
1.05294 + 1.3464 i, -2.60198 i, 0.149446 - 0.448806 i, 

0.149313 

0.0287807 - 2.72544 i, 0.487517 - 1.7952 i, 0.290846 + 1.3464 i, 
0.0287807 - 0.864949 i, -0.619831 i, 0.54 i, 
0.278493 i, -0.277755 i, -0.0641203 i, 0.0641203 i 

0.249771 - 3.14157 i, -3.01053 i, 3.01052 i, 
0.186824 - 2.39256 i, 0.259537 + 2.25 i, 2.3744 i, 
0.186824 - 1.19783 i, 0.54 i, 0.249771 - 0.44882 i, 
0.259537 + 0.442796 i 
0.357565 + 3.13829 i, 0.357656 + 2.69611 i, 0.985138 - 1.7952 i, 
2.15279 i, 0.260208 - 1.7952 i, 1.00394 i, 
-0.98875 i, 0.982836 + 1.3464 i, 0.357565 - 0.445494 i, 
0.357656 - 0.00331156 i 
0.351296 + 2.91719 i, 1.42001 - 1.7952 i, 0.593281 - 1.7952 i, 
0.209089 - 1.7952 i, 0.576676 + 1.3464 i, 0.991553 i, 
-0.988798 i, 0.54 i, 1.41991 + 1.3464 i, 
0.351296 - 0.224399 i 
0.452936 + 2.92003 i, 0.43968 - 2.9168 i, 0.439123 + 2.46302 i, 
- 2.6019 i, 2.15279 i, 0.328707 - 1.7952 i, 
1.04938 i, 0.43968 - 0.673593 i, 0.439123 + 0.229771 i, 
0.452936 - 0.227233 i 
0.625154 - 3.1108 i, 0.624995 + 2.66205 i, -2.6016 i, 
0.830961 - 1.7952 i, 0.238271 - 1.7952 i, 1.69665 i, 
0.825207 + 1.3464 i, 0.54 i, 0.625154 - 0.479587 i, 

0.624995 + 0.0307393 i 
0.755163 + 2.9172 i, -2.60159 i, 1.3168 - 1.7952 i, 
0.550839 - 1.7952 i, 0.199615 - 1.7952 i, 0.530191 + 1.3464 i, 
0.990548 i, 1.31659 + 1.3464 i, 0.54 i, 
0.755163 - 0.224402 i 
1.72415 - 1.7952 i, 0.893009 - 1.7952 i, 0.447923 - 1.7952 i, 
0.176866 - 1.7952 i, 1.70344 i, 0.416589 + 1.3464 i, 
0.891068 + 1.3464 i, -0.988799 i, 0.54 i, 
1.72413 + 1.3464 i 
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E 



Bethe roots Us 



-12.4557 



-9.695 



-8.36086 



-7.49773 



-7.43354 



-6.81246 



-6.75338 



0.484779 - 3.10162 i, 0.411886 + 2.48641 i, 0.106801 - 3.14045 i, 
0.0868008 + 2.51253 i, 0.348418 - 1.88496 i, 0.640811 + 1.25664 i, 
0.106801 - 0.629463 i, 0.0868008 + 0.000739893 i, 0.484779 - 0.668294 i, 

0.411886 + 0.0268653 i 
0.575021 + 2.83026 i, 0.0719595 - 3.14151 i, 0.0613009 + 2.51321 i, 
0.80436 - 1.88496 i, 0.0365392 + 1.72227 i, 0.0365392 + 0.791004 i, 
0.943625 + 1.25664 i, 0.0719595 - 0.628405 i, 0.0613009 + 0.0000685 i, 

0.575021- 0.316986 i 
0.500526 + 2.8322 i, 0.0403235 - 2.64238 i, 0.0383788 + 2.64283 i, 
0.0383947 + 2.38377 i, 0.670337 - 1.88496 i, 0.0403235 - 1.12753 i, 
0.873976 + 1.25664 i, 0.0383788 - 0.129554 i, 0.0383947 + 0.129508 i, 

0.500526 - 0.318927 i 
0.507933 + 2.83108 i, 0.196307 + 3.14091 i, 0.207919 + 2.51414 i, 
0.723661 - 1.88496 i, 0.0062148 + 1.85904 i, 0.899811 + 1.25664 i, 
0.0062148 + 0.654233 i, 0.196307 - 0.627633 i, 0.207919 - 0.000868076 i, 

0.507933 - 0.317806 i 
1.31309 -1.88496 i, 0.0437839 -3.14159 i, 0.0387234 +2.51327 i, 
0.46288 -1.88496 i, 0.00846623 +1.83461 i, 0.635229 +1.25664 i, 
0.00846623 +0.678666 i, 0.0437839 -0.628325 i, 0.0387234, 
1.33352 +1.25664 i 
1.28291 - 1.88496 i, 0.00540715 + 2.56635 i, 0.00749169 - 2.55718 i, 
0.0054066 + 2.46021 i, 0.418583 - 1.88496 i, 0.00749169 - 1.21273 i, 
0.603838 + 1.25664 i, 0.0054066 + 0.0530653 i, 0.00540715 - 0.053073 i, 

1.30642 + 1.25664 i 
0.534826 - 3.13136 i, 0.366787 + 2.83126 i, 0.11739 - 2.52032 i, 
0.0961446 + 2.51326 i, 0.433162 + 1.84874 i, 0.11739 - 1.24959 i, 
0.0961446 + 0.0000111021 i, 0.433162 + 0.664539 i, 0.366787 - 0.317989 i, 

0.534826 - 0.638555 i 



Table 2: The 27 energies and corresponding Bethe roots given by PS?' l \u) for iV = 
3 , s = 1 ,p — 4 , r\ — i7r/5 , cc_ = i] , {3- = 0.35 , Q- = 0.54 , a + = r\ , (3+ = 
0.76,0+ = 0.54 
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E (continued) 



Bcthe roots Uj (continued) 



-6.37707 



-5.9431 



-4.96948 



-4.75362 



-4.29889 



-4.09589 



-4.05775 



-3.93649 



-3.90338 



0.450862 - 3.1307 i, 0.3548 + 2.83363 i, 0.203444 + 3.14135 i, 
0.313088 + 1.93502 i, 0.10853 + 1.81328 i, 0.10853 + 0.699993 i, 
0.203444 - 0.62808 i, 0.313088 + 0.578258 i, 0.450862 - 0.639212 i, 

0.3548 - 0.320357 i 
0.15167 + 3.14151 i, 1.28717- 1.88496 i, 0.152799 + 2.51336 i, 
0.431361 - 1.88496 i, 0.000660592 + 1.88095 i, 0.610477 + 1.25664 i, 
0.000660592 + 0.632323 i, 0.15167 - 0.628234 i, 1.3102 + 1.25664 i, 
0.152799 - 0.0000868625 i 
0.245899 - 3.14157 i, 1.09088 - 1.88496 i, 0.0674913 - 2.51035 i, 
0.056464 + 2.51327 i, 0.236615 + 1.88684 i, 0.0674913 - 1.25957 i, 
1.14613 + 1.25664 i, 0.236615 + 0.626432 i, 0.245899 - 0.628337 i, 

0.056464 

0.118717 + 3.04878 i, 0.118716 - 3.04877 i, 1.08872 - 1.88496 i, 
0.120353 + 1.97697 i, 0.0725214 + 1.8268 i, 0.0725214 + 0.686472 i, 
0.118716 - 0.721137 i, 0.118717 - 0.535507 i, 0.120353 + 0.536301 i, 

1.14446 + 1.25664 i 
0.611823 - 3.10726 i, 0.585252 + 2.48522 i, 0.283735 + 2.82749 i, 
0.545972 - 1.88496 i, 0.00037444 + 1.88259 i, 0.773596 + 1.25664 i, 
0.00037444 + 0.630686 i, 0.283735 - 0.314213 i, 0.611823 - 0.662648 i, 

0.585252 + 0.0280531 i 
0.792496 + 2.51421 i, 0.79911 - 2.5049 i, 0.0126395 + 2.51327 i, 
0.0132493 - 2.51136 i, 0.0132493 - 1.25855 i, 0.0832428 + 1.25664 i, 
0.79911 - 1.26501 i, 0.939738 + 1.25664 i, 0.0126395, 
0.792496 - 0.000938409 i 
0.79092 +2.51422 i, 0.797571 -2.50483 i, 0.0134204 +3.14159 i, 
0.0128183 +1.8868 i, 0.0564763 +1.25664 i, 0.797571 -1.26508 i, 
0.938917 +1.25664 i, 0.0134204 -0.628319 i, 0.0128183 +0.62647 i, 

0.79092 -0.000943414 i 
0.886933 -3.14031 i, 0.0107824 +2.51327 i, 0.0112295 -2.51164 i, 
0.876475 +1.88715 i, 0.725212 -1.88496 i, 0.0112295 -1.25828 i, 
0.074371 +1.25664 i, 0.0107824, 0.876475 +0.626128 i, 
0.886933 -0.629601 i 
0.886042 -3.1403 i, 0.011339 +3.14159 i, 0.875518 +1.88714 i, 
0.723621 -1.88496 i, 0.0108981 +1.88655 i, 0.0516604 +1.25664 i, 
0.011339 -0.628319 i, 0.0108981 +0.626729 i, 0.875518 +0.626131 i, 

0.886042 -0.629609 i 
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E (continued) 



Bethe roots Uj (continued) 



-3.5973 



-2.69459 



-2.39712 



-1.91101 



-1.65102 



-1.3681 



-1.24743 



-0.451049 



-0.278905 



0.336944 -3.14041 i, 0.335901 +2.51212 i, 1.23958 -1.88496 i, 
0.391865 -1.88496 i, 0.000094209 +1.88434 i, 0.568642 +1.25664 i, 
0.000094209 +0.628934 i, 1.26813 +1.25664 i, 0.336944 -0.629498 i, 

0.335901 +0.00115798 i 
0.779303 +2.82839 i, 0.232072 +2.82743 i, 1.05217-1.88496 i, 
0.329315 -1.88496 i, 0.0000460866 +1.88465 i, 0.475658 +1.25664 i, 
0.0000460866 +0.628622 i, 1.10962 +1.25664 i, 0.232072 -0.31416 i, 

0.779303 -0.315113 i 
0.503174 +2.79201 i, 0.504144 -2.78706 i, 0.506179 +2.23281 i, 
0.745784 -1.88496 i, 0.000133753 +1.88409 i, 0.944351 +1.25664 i, 
0.000133753 +0.629187 i, 0.504144 -0.982846 i, 0.503174 -0.278731 i, 

0.506179 +0.280461 i 
0.745717 -3.13841 i, 0.476924 -2.52358 i, 0.476385 +2.51311 i, 
0.712375 +1.88751 i, 0.000108768 +1.88425 i, 0.476924 -1.24633 i, 
0.000108768 +0.629027 i, 0.476385 +0.000165845 i, 0.712375 +0.625767 i, 

0.745717 -0.631505 i 
0.591714 +3.13021 i, 0.594377 +2.52503 i, 1.08954 -1.88496 i, 
0.319176 -1.88496 i, 0.0000183025 +1.88483 i, 0.467049 +1.25664 i, 
0.0000183025 +0.62844 i, 1.14296 +1.25664 i, 0.591714 -0.616939 i, 

0.594377 -0.0117591 i 
1.52581 -1.88496 i, 0.205857 +2.82743 i, 0.710251 -1.88496 i, 
0.25594 -1.88496 i, 1.8849 i, 0.349688 +1.25664 i, 
0.813986 +1.25664 i, 0.628371 i, 0.205857 -0.314159 i, 
1.53319 +1.25664 i 
0.75169 +2.86018 i, 0.738919 -2.84659 i, 0.731307 +2.15882 i, 
0.34774 -1.88496 i, 0.0000213862 +1.88481 i, 0.564622 +1.25664 i, 
0.0000213862 +0.62846 i, 0.738919 -0.92332 i, 0.731307 +0.354451 i, 

0.75169 -0.346911 i 
0.496986 +2.82743 i, 1.4849 -1.88496 i, 0.679662 -1.88496 i, 
0.240941 -1.88496 i, 1.88494 i, 0.326913 +1.25664 i, 
0.789576 +1.25664 i, 0.628339 i, 1.49385 +1.25664 i, 
0.496986 -0.314155 i 
0.901463 -3.10659 i, 0.894513 +2.47987 i, 0.816494 -1.88496 i, 
0.258214 -1.88496 i, 1.88493 i, 0.358314 +1.25664 i, 
0.927674 +1.25664 i, 0.628344 i, 0.901463 -0.663324 i, 
0.894513 +0.0334058 i 
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E (continued) 


Bethe roots Uj (continued) 


0.72331 
1.69087 


1.02774 +2.82771 i, 1.29908 -1.88496 i, 0.601685 -1.88496 i, 
0.212716 -1.88496 i, 1.88495 i, 0.282851 +1.25664 i, 
0.716771 +1.25664 i, 0.628324 i, 1.31854 +1.25664 i, 
1.02774 -0.314435 i 
1.74631 -1.88496 i, 0.943277 -1.88496 i, 0.514308 -1.88496 i, 
0.183505 -1.88496 i, 1.88495 i, 0.238533 +1.25664 i, 
0.614563 +1.25664 i, 0.62832 i, 0.990458 +1.25664 i, 
1.7488 +1.25664 i 
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